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ABSTRACT 

Using the Spitzer Space Telescope, we have obtained rest frame 9-16 /xm spectra of 11 quasars 
and 9 radio galaxies from the 3CRR catalog at redshifts 1.0 < z < 1.4. This complete flux-limited 
178MHz-selected sample is unbiased with respect to orientation and therefore suited to study 
orientation-dependent effects in the most powerful active galactic nuclei (AGN). The mean radio 
galaxy spectrum shows a clear silicate absorption feature (rg./^m — 1.1) whereas the mean quasar 
spectrum shows silicates in emission. The mean radio galaxy spectrum matches a dust-absorbed 
mean quasar spectrum in both shape and overall flux level. The data for individual objects 
conform to these results. The trend of the silicate depth to increase with decreasing core fraction 
of the radio source further supports that for this sample, orientation is the main driver for 
the difference between radio galaxies and quasars, as predicted by AGN unification. However, 
comparing our high-z sample with lower redshift 3CRR objects reveals that the absorption of 
the high-z radio galaxy MIR continuum is lower than expected from a scaled up version of lower 
luminosity sources, and we discuss some effects that may explain these trends. 

Subject headings: galaxies: active; galaxies: high-redshift; infrared: galaxies 



1. Introduction 

There is strong evidence that many lumi- 
nous active galactic nuclei (AGN) are surrounded 
by dust in a torus-like geometry leading to 
orientat ion-dependent obscuration of the nuclear 
region (lAntonucci fc Miller I Il985l: iBarthell [l989t 



Antonuccilll993l : lUrrv &: Padovanilll995l ). Despite 

the general consensus on the AGN unification 
principle, it is still a matter of debate whether 
obscuration is indeed the dominating effect on the 
observed differences between radio galaxies and 
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quasars. The question remains if every obscured 
(type-2) AGN has a (type-1) AGN hidden by a 
"torus" or if some of the observed differences may 
arise from other effects, for instance different evo- 
lutionary states, different star forming activity in 
the hosts, or different accretion rates. 

Orientation-dependent effects can only be 
tested with type-1 and type-2 AGN samples 
matched in isotropic emission. Low-frequency 
(meter wavelength) radio-selected AGN are partic- 
ularly attractive for studying orientation-dependent 
properties at other wavelengths because the inte- 
grated emission from the radio lobes is optically 
thin and essentially isotropic. This paper there- 
fore focuses on powerful double-lobed radio galax- 
ies (type-2) and steep-spectrum qu asars (type- 



1) ta ken from the 3CRR catalogue (jLaing et al 
19831). a subset of the de eper and larger 3GR sam- 



ple (jSpinrad et al.lll985l ). The 3CRR sample pro- 
vides a complete, flux-limited sample of radio-loud 
objects selected at 178 MHz. Mid infrared (MIR) 
and far infrared (FIR) observations of orientation- 
unbiased samples are valuable tools for the study 



of unification schemes. They probe the intrinsic 
energy budget of the nuclear regions where the 
optical/UV radiation from the accretion disk is 
absorbed by the surrounding dust, which then re- 
emits the absorbed energy at infrared wavelengths. 
Thus, the infrared might serve as a "calorimeter" 
for the radiative output in the nucleus. How- 
ever, while the infrared radiation has the ability 
to penetrate the obscuring dust, it will still suffer 
some (wavelength-dependent) extinction, which 
may yield important clues as to the nature of the 
obscurer itself. 

At lower redshift (z < 0.5) the large number 
of radio galaxies compared to quasars is incon- 
sistent w ith orientati on-dependent models (e.g.. 
Fig. 1 of ISingallll993l) . It has been argued that 
this might be due to an additional population of 
low-luminosity radio galaxies which falls below 
the selection limit in typical higher-redshift, low- 
frequency-selected samples. These low-luminosity 
radio galaxies often show quite different multi 
wavelength properties with moderately absorbed 
(and less luminous) accretion-rel ated X-ray emis- 
sion (e.g. iHardcastle et al.l 120091) . optical sp ectra 
of low ionization Ce.g. lButtiglione et al.ll2009l). and 
weak MI R emission fMeisenheimer et al.l 12001 : 
Ogle et al. 2006; van der Wolk et al. 201^ 

At intermediate redshifts (0.5 < z < 1.0) 
numerous studies — often utilizing IR data — 
show that the unification picture works well for 
radio galaxies a nd quasars (e.g. iBarthell Il989t 



Haas et al.l 120041 ). To date, a large number of 



z < 1 SCR objects have be en studied in the 
MIR and FIR with ISO (e^ Ivan Bemmel et al 



2000llMeisenheimer et al 2001t and as compiled 
by Siebenmorgen et al.l 12004 and by Haas et al 



2004 ) and with 6'iPzteer ( e.g.lShi et al.H2005|lHaas et al 
20051 : iQgle et al.ll200fil : ICL 



earvet al.ll2007l ). While 



such studies were vital to establish the unifica- 
tion picture at 0.5 < z < 1.0, some of this work 
also indicated that on average quasars seem to 
be a factor of a few more luminous in the infrared 
than radio galaxies when normalized by their total 
radio luminosity. To explain this effect, it was ar- 
gued that non-thermal contributions of a beamed 
radio core increase the infrared fiuxes of quasars in 
addition to the wavelength-depen dent obscuration 
of the torus in radi o g alaxies (e.g 



Hes et all 199 



Clearv et al.l l2007l). JMeisenheimer et al.l pOOll 



ence between quasars and radio galaxies but did 
not model it in detail (though the simplest expla- 
nati on — obscuration -- was b riefly mentioned), 
and ISiebenmorgen et al.l (|2004l) have interpreted 
the difference in accordance with an obscuration 
scenario. 

At high redshifts (z > 1.0) the situation is less 
clear. For the SCR sources at these redshifts, sen- 
sitive FIR (60-200 /im) observations have yet to be 
made, but S.6-24/im photometry of the complete 
high-z part (1.0 < z < 2.5) of the SCR sample 
has been obtained with Spitzer. The rest frame 
2-10 /im spectral energy distributions (SEDs) pro- 
vide evidence for dust emission which is a factor 
3-5 w eaker in type- 2 A GN compared to type-1 
AGN ( Haas et al.l 12008 ). The mean radio galaxy 



SED is consistent with a reddened version of the 
mean quasar SED, leaving little room for a signifl- 
cant contribution of beamed non-thermal emission 
to the infrared SEDs of quasars, at least for that 
sample. 

Apart from the unification issues, it is vital to 
further explore the largely unknown MIR proper- 
ties of high redshift radio-loud AGN and to ex- 
pand the luminosity range and potentially the 
AGN-host contrast studied so far. Therefore, 
we here analyze Spitzer spectra of the quasars 
and radio galaxies in the SCRR sample at red- 
shift 1 < z < 1.4. We use a ACDM cosmology 
with Ho = 71 km s'^ Mpc-\ fl^ = 0.27, and 
riA = 0.7S. 

2. Data 



Our sample comprises all objects of the SCRR 
catalogs with 1 < z < 1.4j. This selection includes 
20 sources in total, 11 quasars and 9 radio galaxies. 
Some basic properties of the sources are presented 
in Table [TJ The redshift range 1 < z < 1.4 was 
selected to cover the rest frame wavelengths 9- 
16 /im when using the LLl sp ectral module of the 
Infrared Spectrograph (IRS, JHouck et al] |2004|) 



on b oard Spitzer Space Telescope ([Werner et al 
20041 ). For 1 < z < 1.4, the spectral win- 



and iHaas et al.l (|2004l ) saw the MIR/FIR differ- 



^ http://3crr.extragalactic.info/ 

^Due to an updated redshift llGrimes et a l.'l2005'). 3C325 
is included. 4C16.49 (z = 1.296: lAarsTT al. 2005) should 
in principle be part of the sample, but at the time the 
observations were defined the NED redshift was z = 1.88, 
and therefore no data were obtained for this source. 



Table 1: The Sample. 



Object 


z 


-Pi 78 MHz 
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pcorc 
-^5 GHz 
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i^l5 


/^m 


Tg. 7^111 






(Jy) 




(mJy) 




(mJy) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 








Q 


uasars 










3C 068.1 


1.238 


14.0 


0.80 


1.1 


4 


9.3 


0.7 


-0.05 


3C181 


1.382 


15.8 


1.00 


9.0 


8 


8.7 


0.5 


0.19 


3C186 


1.063 


15.4 


1.15 


15.0 


7 


8.2 


0.6 


-0.12 


3C190 


1.197 


16.4 


0.93 


73.0 


7 


13.4 


1.5 


0.60 


3C204 


1.112 


11.4 


1.08 


26.9 


4 


8.3 


0.5 


-0.13 


3C208 


1.109 


18.3 


0.96 


51.0 


4 


5.8 


0.5 


-0.27 


3C212 


1.049 


16.5 


0.92 


186.0 


1 


15.5 


0.9 


-0.11 


3C245 


1.029 


15.7 


0.78 


1219.0 


12 


28.3 


1.9 


-0.02 


3C 268.4 


1.400 


11.2 


0.80 


45.0 


11 


15.6 


0.8 


0.05 


3C287 


1.055 


17.8 


0.42 


3280.0 


9 


7.1 


0.6 


-0.23 


3C325 


1.135 


17.0 


0.70 


2.4 


5 


4.5 


0.3 


-0.21 


Radio Galaxies 


3C013 


1.351 


13.1 


0.93 


0.18 


2 


5.6 


0.9 


0.85 


3C065 


1.176 


16.6 


0.75 


0.52 


2 


2.2 


0.2 


0.30 


3C252 


1.105 


12.0 


1.03 


2.2 


5 


11.8 


1.0 


0.61 


3C266 


1.272 


12.1 


1.01 


0.4 


6 


2.9 


0.4 


0.72 


3C267 


1.144 


15.9 


0.93 


1.87 


10 


6.0 


0.5 


0.77 


3C324 


1.206 


17.2 


0.90 


0.17 


3 


8.3 


0.7 


1.71 


3C356 


1.079 


12.3 


1.02 


1.1 


5 


7.4 


0.6 


0.91 


3C368 


1.132 


15.0 


1.24 


0.22 


3 


8.1 


0.7 


1.74 


3C 469.1 


1.336 


12.1 


0.96 


2.4 


10 


8.3 


1.3 


1.74 





















Note. — (1) Object name; (2) Redshift; (3) Observed 178 MHz flux in Jy; (4) Spectral index (5^ oc i/~") between (observed) 
178 and 750 MHz as taken from the 3CRR web page; (5) Flux density of t he radio core in m Jy; (6) References for the radio-core 
flux densities. They refer to observed frame 5 GHz except for data from lBest et al.l l ll997t) . which were taken at 8.4 GHz; (7) 
Flux density measured at rest-frame 15 ^m in mJy; (8) Icr uncertainties on the 15 /xm flux density; (9) Screen opacity relative 
to the mean quasar spectrum as determined in i|3.2l 



References. — (plAkuj or ct al. 1991: (2) Best ct al. 1997; (3) B est et allll998l: f4V Bridlc ct al."l994': (5) 'Fcrnini ct alj ll997l; 
fel iLiu et al.lll992l : ( 7) iLiidkc ct al. 1998 ; (8) Mantovani et al. 19921: (91 [Pearson et al..il985; (10) Pedelty ct al. 1983 ; (11) 
iReid et al.lll995l : (12') ISaikia et al.lll990l 



dow of the IRS/LLl (19.5 — 36.5 /im) includes 
redshifted low- and high-excitation atomic emis- 
sion lines [Ne II]A12.8^m, [Ne V]A14.3^m, and 
[Ne III] A15.5/im. It also covers potential 11.3 /im 
PAH emission from star formation and the silicate 
feature at 9.7 /xm. 

18 objects were observed for a total of 2880 
seconds each with IRS in the LLl module (PID- 
40314, PI M. Haas). The remaining objects, 
3C 356 (PID 03349, PI R. Siebenmorgen) as well as 
3C325 (PID 00074, PI J. Houck), were observed 



previously in LLl for a total of 960 s each. The 
individual frames were averaged for each nod po- 
sition separately, and the resulting images were 
used to subtract the background from the other 
nod. The final background-subtracted images 
were cleaned of residual rogue pixels and cosmic 
rays using irsclean. The SPICE software was used 
to calibrate and extract one dimensional spectra. 
The extracted spectra from the two nod positions 
were then averaged into a single spectrum. 

The comparison of our spectroscopy with pho- 
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Fig. 1. — Spectra and photometry for the mdividual objects in the rest frame of the sources in arbitrary 
linear units of lyL^,. SEDs have been offset verticaUy for clarity, but 15 fim flux densities can be found in 
Table 1. Quasars are shown on the left and radio galaxies on the right. The sources are ordered according 
to their 15 fim luminosities with i^£i/,i5^m increasing from bottom to the top. The filled symbols show the 
observed photometry while the open symbols represent the interpolated photometry used for constructing 
the average SEDs (Fig. [2]). To the right of each SED we give the object name, the redshift, and the rest 
frame luminosity J^Li/.is^m in ergs~^. 
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Fig. 2. — Average 9 — 16 /im spectra of the high-z 3CRR sources supplemented by 1 — 8/im photometry. 
The solid lines represent the average spectra while the symbols indicate the average photometry. The SEDs 
are normalized to the isotropic low-frequency radio luminosity. The black dashed line represents the average 
quasar spectrum reddened by a dust screen of Ay — 20 mag without any additional scaling. The shaded 
areas indicate the la dispersion. 



tome tric measurements in the MIPS (|Rieke et al. 
20041 ) 24 /im filter reveals a good agreement ex- 
cept for 3C13, 3C65, 3C68.1, and 3C267, all of 
which show a flux deficit in the spectra. Checking 
the pointing for the observation reveals that for 
3C 13, the spectroscopic slit fell slightly off-center 
thus missing part of the flux. We applied a multi- 
plicative factor of 1.8 to the spectrum to bring it 
to the level of the photometric measurement. For 
the radio galaxy 3C 267, the MIPS 24 ^m image 
reveals a faint source to the north-west which is 
also partly covered by the slit. The distance of 
this faint source from the radio galaxy is close to 
the nod throw along the slit. While this config- 
uration can potentially lead to an overcorrection 
when the two nod positions are subtracted from 
each other to remove the background, we cannot 
identify the faint source on averaged 2-dimensional 
spectral frames at the individual nod positions. A 
multiplicative factor of 1.2 is needed to match the 
spectrum with the photometry. For the other two 
sources, 3C 65 and 3C 68.1, no immediate cause for 
the fiux difference could be identified. For these 



objects we applied correction factors of 1.5 and 
1.3, respectively. These corrections do not affect 
the basic conclusions drawn here, but we mention 
that wavelength-dependent slit losses could result 
in overestimating the silicate depth in 3C 13 (sec- 
tion [Ol. 



2.1. Creating Average Spectra and SEDs 

All objects have been securely detected, and 
continuum features such as silicates in emission or 
absorption can be observed in individual sources 
(see Table [T] for measured rest-frame MIR fiuxes 
and FigH] for the observed spectra) . Because the 
quality of the spectra is often limited, we created 
separate average spectra for the radio galaxies and 
quasars. First, the individual spectra were shifted 
into the rest frame of the source, interpolated onto 
a common wavelength grid, and then the inter- 
polated spectra were transformed into vL^. To 
avoid the average spectra being dominated by the 
most MIR luminous objects, the spectra were nor- 
malized in i^ii/,i5pm before averaging. While this 



normalization is well suited to emphasize differ- 
ences in the spectral shape, information about the 
absolute MIR luminosity is removed from the av- 
erage spectra. This information, however, may be 
recovered using the extended radio luminosity of 
the sources. 

Due to their selection, the quasars and radio 
galaxies in this sample span the same range of 
radio luminosities. The actual rest frame mean 
values are (i^Li^^iysmhz) = (1.9 ±0.6) x 10''^erg/s 
for the quasars and {i^Li,^nsMHz) = (2.0 ± 0.4) x 
10^"* erg/s for the radio galaxies. First we used the 
15 /im normalization factors to scale the individual 
^^ I/, 178 MHz values. Then we normalized the aver- 
age IR SEDs of each class to have the same aver- 
age lyL^.nsMHz of j^L^,i78mhz = 2.0 x 10*^ erg/s, 
which is the mean rest frame 178 MHz radio lumi- 
nosity of the complete sample. 

In addition to the spectra, observe d photome- 



try a t 3.6, 4.5, 5.8, 8.0, 16, and 24 /im (jHaas et al 



I2OO8I ) was used to construct average near-infrared 
(NIR)/MIR SEDs for the sample. The photo- 
metric measurements were shifted in wavelength 
into the rest frame of the individual sources and 
then spline interpolated onto a common wave- 
length grid. Figure [1] shows the observed as well 
as the interpolated photometry together with the 
spectra. After transforming the interpolated val- 
ues into vL^ we applied the same scaling factors 
as for the spectra and then averaged the photom- 
etry. After normalization, all the quasar SEDs 
are very similar to each other in shape (except 
for the strength of a NIR bump, see section [3?T1) . 
while the radio galaxy SEDs tend to be more di- 
verse (Fig. [ p. Th is has already been explored by 
Haas et al.l (|2008[) for their larger and fainter sam- 



ple (which includes the current sample), and the 
trends are similar. Despite the diversity in the ra- 
dio galaxy SEDs, we did not exclude any object 
from the averaging process. 

3. Results and Discussion 
3.1. Average SEDs 

Figure [2] shows the average spectra and pho- 
tometry for the radio galaxies and quasars. The 
combination of the deep silicate absorption at 
9.7 /zm and the overall lower flux levels indicate 
that the radio galaxies on average suffer consid- 
erable absorption throughout the entire MIR con- 



tinuum relative to quasars. Applying a redden- 
ing screen of Ay = 20 mag to the average quasar 
spectrum yields a very good match with the radio 
galaxy spectrum in the continuum level as well as 
in the shape of the silicate absorption featured 
From the MIR point of view, this strongly sup- 
ports the idea that the quasars and radio galaxies 
in our sample are in fact intrinsically identical with 
the former being essentially unobscured while the 
latter suffer from extinction. The good match of 
the average radio galaxy spectrum to the reddened 
quasar spectrum also implies that the intrinsic ra- 
dio galaxy spectrum is not flat. In particular, the 
average radio galaxy must have silicate emission 
underlying the absorption feature (see section !??^ . 
We here assumed a screen scenario where the 
absorbing material is placed in front of the emit- 
ter. However, in the clumpy torus models which 
in recent years were able to successfully explain 
the dust emission in local Seyfcrt galaxi es (e.g 
H5nig et al.l l2006t lEhtzuiJ [2OO8; Ncnkova et al 



20081 ) ■ the MIR-emitting zone is not located at 



a single radius but is spread over a range of radii, 
thus essentially creating a mixed-case scenario for 
parts of the absorber. In general, screen extinc- 
tion underestimates column density by different 
amounts at different wavelengths. In addition, the 
bulk of the MIR emission will originate at mod- 
erate depth within the obscuring structure thus 
suffering (far) less extinction than radiation from 
the accretion disk or from hot dust located near 
the accretion disk. Thus, depending on the com- 
plexity of the geometry, the mean Ay determined 
from the spectroscopy represents a lower limit and 
will in most cases underestimate the total column 



densi ty by a large factor (see e.g. iLevenson et al 



20071) 



A notable feature in the average quasar SED 
is the prominent NIR bump at '^2-5 /xm (Fig. [2]). 
This component is, to different degrees, present 
in the SEDs of essentially all individual quasars 
(Fig. [T]). Such NIR bumps have been previ- 
ously rep o rted f or individual type-1 AGN (e.g 



'Barvami; 

;2006 : Riffel et al 



i^ [l98i 



Rodriguez- Ardila fc Mazzalay 
20091 ) and can be seen also 



in larger samples of quasars and in composite 



^We here use the MIR extinction curve from 
IChiar fc Tieleng ||2006| ) normalized to the peak of 
the sihcate feature at 9.7/im and Ay/T9.7fini = 18. 



SEDs fe.g. lEdelson fc Malkaij llQSd: lElvis et al 



19941 : iGallagher et alJ 120071 ) . They are generally 



identified with graphite dust near the sublima- 
tion temperatur e located very close to the ac - 
tive nucleus fe.g. lBarvainislll987tlMor et al.ll2009[) . 
There is some indication that the strength of this 
feature increases with i ncreasing visible/UV lu- 
minosity of the quas ar ( Edelson fc Malkanlll986t 
[Gallaghe 7eraDl2007t ) . 

In contrast to the quasars, the radio galaxies 
show no noticeable 2~5 /im bump. Screen redden- 
ing the average quasar SED with Ay = 20 mag 
leaves an excess for A < 8/im and especially at 
the location of the NIR bump (FiglH top). This 
feature is either missing in the radio galaxies, or 
it might suffer higher extinction than the longer- 
wavelength continuum. 

(|2n08h 



Haas et al. 



For their larger sample, 
found that radio galaxy SEDs between 2 and 
10 /zm are on average consistent with reddened 
quasars, similar to our result from the spec- 
troscopy at 10-16 /im. Because their sample of 
SCR objects contains a large number of fainter 
and higher-redshift sources compared to the sam- 
ple of 3CRR objects studied here, we performed 
the same analysis as described by Haas et al. but 
limiting the objects to our 20 sources. Using their 
approach, it would require a screen extinction 
of Ay ^ 40 mag to match the average 2-10 /im 
SEDs of radio galaxies and quasars in our sample. 
This value is higher than the value determined 
from the MIR spectroscopy alone {Ay ^ 20 mag, 
Tg.7i_im ~ 1.1 mag) as discussed above. 

For the discrepancy at A < 8 /im, several pos- 
sibilities can be imagined. One example is in- 
accuracies in the MIR extinction curves. Dif- 
ferent extinction curves show differences in the 
continuum at ^8 ^m and in the silicate fea- 



ture ( e.g. iMathis et al 



19851 : IChiar fc Tieleii? 



1983t iRieke fc Lebofskv 



20061 ). 



and an incorrect 
ratio between continuum and silicate extinc- 
tion could be part of the problem. However, 
considering the good fit in the silicate depth 
and in the continuum extinction for the spec- 
tra (Fig. [H tojo), the e xtinction curve used here 



(jChiar fc Tielens 
longer wavelengths, 
tinction curves (e.g. , 



20061 ) works well at least at 

Moreover, publish e d ex- 

Indebetouw et al.l l2005t 



Nishivama et al.ll2009l ) are all similar and equally 
featureless shortwards of the silicate feature. With 



30 f 



■So 20 - 
U - • • 



^ 10 
> 



oi 

30 f 




'^BB,').7^lll — O'i 

T„„=1324K 



20 



^ 10 
> 




ol 

30 [ 



T„„=I435K 



20- 



• • 



^ 10 
> 




5 10 15 

rest wavelength (|J,m) 



20 



Fig. 3. — Average quasar and radio galaxy SEDs 
compared to dust screen reddening (dashed lines) . 
No additional scaling was applied. The data are 
those shown in Fig.[2]but with uncertainties omit- 
ted for clarity. Top: Quasar SED reddened us- 
ing a single rg.y^m value. Middle: Quasar SED 
reddened with two components, one affecting all 
wavelengths and the other, larger one affecting 
only the single blackbody component represent- 
ing the 3-5 /im bump. Bottom: Quasar SED with 
the blackbody component representing the 3-5 /im 
bump subtracted and the remainder reddened by 
a single screen value. 



curves such as these, it is not possible to match 
the NIR and MIR simultaneously in a screen con- 
figuration affecting all wavelengths equally. 

An alternative solution can be found by includ- 
ing in the fits a separate blackbody component 
representing the NIR bump. Although a solid the- 
oretical foundation for the origin and the prop- 
erties of the NIR bump is still missing, it can 
be well d escribed as a single temperature black- 
body fe.g. lBarvaiifiillQSTtlMor et al.ll2009t ). If the 
bump comes from dust, the necessary high tem- 
perature places this dust component close to the 
nucleus, and the radiation from this component 
would have to penetrate a larger dust column than 
the rest of the continuum emission. Assuming a 
black body component from nuclear hot dust suf- 
fers more extinction than the rest of the IR contin- 
uum, screen-reddening of the average quasar SED 
is consistent with the average radio galaxy SED 
(Fig. [21 middle), and the derived temperature of 
the NIR bump blackbody (T ~ 1300 K) is in the 
range expe cted from previo us studies in the liter- 
ature (e.g. iMor et al.l 120091 ). The major require- 
ment in this scenario is that the screen opacity 
for the NIR bump component (Tg.y^ni « 6.1) must 
be significantly higher than the opacity affecting 
other components (rg.y^m ~ 0.9). 

It could also be suggested that the NIR bump 
is absent in most radio galaxies. Subtracting a 
single-temperature blackbody from the average 
quasar spectrum before screen reddening the re- 
maining continuum gives a good fit as well (Fig. [3l 
bottom). However, this scenario would imply 
that there are intrinsic differences between ra- 
dio galaxies and quasars, where the latter host 
an additional dust component close to the nu- 
cleus. Such a scenario seems unlikely for the ma- 
jority of the objects. Visible-light imaging and 
polarization studies show that at 0.7 < z < 2.5, 
essentially all luminous ra dio galaxies show the 
so-called alignment ef fect ([Chambers et al.lll987 : 
iMcCarthv et al.l [l987 ) , and much of this emis- 
sion turns out to be scattered l i g ht from a 



hidden quasar (e .g. [Cimatti et al.l [1993L [1997 : 



Vernet et al.ll200l[ ). Thus, at higher redshifts most 



if not all radio galaxies host a hidden quasar, ar- 
guing for the similarity of both types of objects. 
This also strongly argues that the main difference 
between these two types of objects is in the orien- 
tation of their nuclear structures with respect to 



us. Further clues on the nature of the NIR bump 
might be provided by detailed radiative transfer 
modeling of the SEDs of the sources (F. Heymann, 
in preparation) . Regardless of the reason why the 
NIR bump is not seen in radio galaxies, the esti- 
mated value for Tg.7^„i changes very little. That 
is because it is determined mostly by the (rest 
frame) 10-13 /im slope and secondarily by the 
(rest frame) 8/im photometry and hardly at all 
by shorter-wavelength assumptions or data. The 
rest of this paper uses only r9.7^ni, and the NIR 
difference does not affect our main conclusions. 

3.2. Individual Spectra 

3.2.1. Fitting the Silicate Feature 

A striking feature in Figure[2]is the deep 9.7 /im 
silicate absorption observed for the radio galaxies. 
The long wavelength side of this feature is well 
measured in the spectra, and the shorter wave- 
length photometry data confirm the presence and 
depth of this feature. The silicates alone provide 
clear evidence that in the mean radio galaxy, a 
substantial amount of absorbing dust lies between 
the MIR continuum emitter and the observer. 

For further analysis we determined the depth 
of the silicate feature in the individual objects 
by fitting the observed spectrum and photome- 
try with an input spectrum reddened by a dust 
screen. Considering the uncertainties in treating 
the NIR bump component, we limited the fits for 
the individual objects to A > 7 /z m. We again used 
the ex tinction curve presented bv [Chiar fc Tielens 
( 20061 ) and assumed that the intrinsic continuum 
of every source resembles the average quasar SED, 
including silicates in emission. In order to param- 
eterize the average quasar SED, it was fit by a 
linear combination of three blackbodies of vari- 
able temperature. One of these blackbodies was 
multiplied by the MIR extinction curve to include 
a silicate emission feature We do not assign any 
physical interpretation to this fit as it was only 
used to yield a smooth representation of the av- 
erage quasar continuum. This smooth continuum 
was then taken as the intrinsic continuum for the 
reddening fits to the individual objects. 

Figure [4| shows two examples for our fitting 
procedure, and the derived values for rg.y^m for 
all objects are given in Table [1] As in the ra- 
dio galaxy average spectrum, all individual galaxy 
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Fig. 4. — Two examples illustrating the fitting of 
the silicate feature for a quasar (top) and a radio 
galaxy {bottom). The solid line and filled circles 
show the spectroscopy and interpolated photom- 
etry, respectively. The dashed line shows the fit 
which is a scaled and screen reddened version of 
the average quasar SED. 



data longwards of ^3 fim can be fitted by intro- 
ducing stronger reddening of a NIR bump (or as- 
suming such a component is absent). In the ma- 
jority of cases, variations of the screen rg.r^m with 
other fit parameters are smaller than ±0.1. There- 
fore we assign this value as an uncertainty to the 
individual rgy^m screen values. 

Silicates appear in absorption in all radio galax- 
ies and, with the exception of 3C 190, in emis- 
sion in all quasars. Thus the average spectra 
accurately represent the two populations. Un- 
like Galactic sources, however, the silicate emis- 
sion feature in the average quasar spectrum peaks 
at ^11 fiui. Such a shift in the peak wavelength 
of the silic ate emission is freq uently obser ved in 



AGNs fe.g.lSiebenmorgen et ah 2005: Stur m et aL. 
200,4 iHao et al.ll20n5l : IShi et al.ll20n6h . Under the 



assumption that the emission is optically thin, 
the shift can be explained as a result of fold- 
ing the silicate opacity curve with the (steeply 
rising) black body curve of the silicate grains 
( Siebenmorgen et al.l 120051. Lower temperatures 
would then shift the peak of the silicate feature 



to longer wavelengths. Alternative explanations 
for the shift of the silicate emission peak wave- 
length in AGN include variations in the d ust grain 
size distribution (e.g. ISturm et al. 20051). effects 



from dust chemistry fe.g. lMarkwick-Kemper et al 



2007t Li et al. 20081). or radiative transfer effects 
fe.g. lNikutta et al.ll2009l ). 

The only quasar not showing silicate emission 
is 3C 190, which shows a shallow silicate absorp- 
tion feature instead (Table 1). 3C 190 also has 
the reddest MIR spectrum of all quasars stud- 
ied in this paper. Notably, 3C 190 is a compact 
steep spectrum (CSS: projected linear radio sizes 
smaller than 15kpc) source in the radio, a class of 
objects which may represent young or frustrated 
sources. This class may show different MIR prop- 
erties compared to large double-lobed sources (P. 
Ogle, in preparation). However, the two other 
CSS sources in our sample (the quasars 3C 186 and 
3C 287) appear very similar to the other quasars 
in their MIR properties, including obvious silicate 
emission. 

The radio galaxy 3C65 stands out by having 
a very small T9.7p,m value compared to the other 
radio galaxies. It was also one of the sources where 
we applied a multiplicative factor to match the 
spectrum with the higher photometric 24 /xm flux. 
As discussed in §2, pointing errors are more likely 
to cause an overestimate than an underestimate of 

Tg.T/iin- 

3.2.2. Silicate feature and orientation 

The observed extinctions can be connected with 
the orientations of the sources by analyzing the 
radio properties of the objects (Table [T]). Due 
to relativistic beaming effects, the 5 GHz flux of 
the compact radio core in our objects depends to 
some extent on the orientation of the radio jet. 
Therefore the ratio of the radio core fiux and the 
(isotropic) extended fiux can be used as an indi- 
cator of the approximate inclinati on of the radio 
source (e.g.. lOrr fc Brownd ll982J ). Because the 
radio jet is thought to be oriented along the po- 
lar direction with respect to the obscuring struc- 
ture, one might expect an anti-correlation between 
the core fraction and the extinction in the op- 
tical and MIRo This has been demonstrated to 



*We here determine the core fraction as the ratio between 
the 5 GHz radio core flux (Tab. [l} and the 5 GHz total flux 
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Fig. 5. — Radio core fraction plotted versus the 
depth of the siUcate feature. The core fraction 
was determined as the ratio between the 5 GHz 
core flux and the 5 GHz total flux. 

be the case from the optical spectra of quasars 
which show redder continua (and larger Balmer 
line ratios) with decre asing core fraction (e.g. 
Baker fc Hunsteadlll995() . Figure [5] shows the core 
fraction of our objects as a function of rg.y^m. All 
sources in this figure with high silicate extinction 
have weak cores, while most sources with silicate 
emission show stronger cores. This supports an 
orientation-based explanation for the differences 
between the MIR properties of radio galaxies and 
quasars in the current sample. 

3.2.3. Silicate extinction and 15 fim luminosity 
While earlier studies (e.g. iHaas et al.l 120081 ) 



showed that the MIR difference between quasars 
and radio galaxies is consistent with absorption, 
the actual amount of absorption was still un- 
known. This necessary piece of information can be 
derived from the silicate feature depth. The aver- 
age spectra in Figure [2] already indicate that radio 
galaxies are on average less MIR luminous than 
quasars and that this is consistent with absorp- 
tion. Figure |6] shows the comparison for individ- 
ual objects. Quasars and radio galaxies show con- 
siderable overlap in their 15 /im luminosities, but 
the depth of the silicate feature separates the two 
classes almost completely. Radio galaxies occupy 
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Fig. 6. — Low-frequency normalized luminosity at 
15 /im plotted versus the depth of the silicate fea- 
ture. The arrow shows the effect of a dust screen 
with A y = 20 mag, assuming the IChiar fc TielensI 
( 20061 ) extinction curve. Dereddening the radio 
galaxies along the direction of the Ay arrow shifts 
them to the same normalized Li5^m-range as pop- 
ulated by the quasars. 
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taken from lLaing fc Peacockl | |198C| ). We used the observed 
fluxes and did not apply a K-correction. 



a larger range of observed Tg.y^m than quasars but 
with systematically higher values. The striking 
feature of Figure [6] is that dereddening will match 
up the radio galaxies with the quasars in both the 
silicate strength and radio-normalized MIR lumi- 
nosity. 

The characterization of radio galaxies as red- 
dened quasars appears to be in conflict with re- 
sults at lowe r redshifts . For o bjects with strong 
radio cores, iHes et al.l (|l995[ ) and ICleary et al.l 
(|2007l ) have proposed that beamed synchrotron 
emission can make sig nificant contr i bution s even 
at MIR wavelengths. IClearv et al] (|2007t ) com- 
pared the MIR continuum properties of a sample 
of 0.5 < z < 1.0 radio galaxies and quasars also 
selected from the 3CRR catalogue. They found 
that, normalized to their 178 MHz luminosities, 
quasars are on average ~4 times more luminous 
than radio galaxies at (rest) 15 /itm. Based on SED 
fitting with a spherically symmetric dust model, 
Clearv et al.l assigned half of this difference to non- 
thermal emission from a beamed radio core in the 
quasars while the remaining half was thought to 
arise from MIR absorption in the radio galaxies. 
This combination agrees with standard AGN uni- 
fication schemes. 
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Except for cases borderline to flat spectrum 
quasars such as 3C245, our steep-spectrum ob- 
jects lack such strong cores, and enhanced MIR 
emission due to non-thermal radiation seems un- 
likely. If the MIR continuum flux of the quasars 
in our sample were enhanced by synchrotron emis- 
sion, they would lie systematically at higher lu- 
minosities than the de-reddened radio galaxies 
(Fig. [6]). Furthermore, we would not expect the 
reddened quasars to fit the radio galaxies in both 
continuum level and silicate feature strength with- 
out any corrections for non-thermal contributions 
(Fig. [2]). Thus, while a non-thermal component at 
some level is not entirely ruled out, the results for 
our sample do not require one. 

3.3. Radio- normalized 15 /im luminosity 

Our data set expands the comparison between 
low and intermediate redshift (and luminosity) ra- 
dio sources to more powerful sources. A compari- 
son provides new insight on the luminosity depen- 
dence of the radio- normalized 15 /im dust emis- 
sion. 

3.3.1. Comparison data 

We calculated i?dr = vL^^i^j ^ai/ vLu,n?,MYiz for 
our objects as well as for two comparison samples 
at 0.5 < z < 1.0 and at z < 0.5. For the compar- 
ison samples, we used all FR-II (jFanaroff fc Rilev 
19741 ) radio sources in the 3CRR catalog observed 
spectroscopically with Spitzer. Spectroscopic MIR 
data are available for 42 out of 49 sources (86%) at 
0.5 < z < 1.0 and for 40 out o f 60 source s (67%) 
at z < 0.5. We used data from lOgle et al.l (12006.) 
Clearv et all (|2007t) . and iHardcastle et al.l (|2009l ) 
or retrieved the spectra directly from the Spitzer 
archive and, after standard data reduction, deter- 
mined the MIR fluxes. Rest frame radio and MIR 
luminosities were calculated in the same manner as 
for our high-redshift sample using observed spec- 
tral indices from the 3CRR webpage. Figure [7] 
shows i?cir versus redshift. 

Especially at lower redshift, a population of 
MIR-underluminous radio galaxi es has been found 

200 it IShi et all [20o5 



Meisenheimer et al. 



Ogle et al.M2006[ ). Such objects also have quite 
different multi- wavelength proper ties compared to 
MIR-luminous radio galaxies (e.g. lButtiglione et al 
2009t IHardcastle et al.l 120091 ) . and they may rep- 
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Fig. 7. — log i?dr plotted versus the redshift z. i?dr 
is defined as i^Lj/,i5pm/'^-^i/,i78MHz, and 15 //m and 
178 MHz refer to rest wavelength and frequency. 



resent a population of radio galaxies with low (ra- 
diative) accretion power. We therefore present the 
statistics below excluding all the "MIR -weak" ob- 
jects with i^Li/,i5/ini < 8.0 X W^^erg/s (jOgle et al. 
20061 ). This excludes one source at 0.5 < z < 1.0 
and 17 sources at z < 0.5. All three samples at 
the different redshifts include a small number of 
CSS sources. Excluding those sources would not 
affect the statistics in any significant way. 



3.3.2. Luminosity dependence 

Figure [5] shows R^^ as a function of radio lumi- 
nosity -/jiTsMHz- The most striking features of this 
figure are: 

1) The i?(ii. of quasars decreases steadily with 
increasing Lirg mhz- This decline is seen over 
almost three orders of magnitude in Li78 mhz 
and over the full redshift range (Fig. [7]) cov- 
ered here. 

2) For the radio galaxies (excluding the MIR- 
weak sources), i?dr at low to intermediate 
-^178 MHz decreases similar to the quasars. 
However, for Li78MHz ^ 2 x 10'^^ erg/s/Hz, 
i?dr flattens or even increases towards the 
highest luminosities. This leads to an ap- 
parent i?dr convergence of quasars and radio 
galaxies in the highest luminosity/redshift 
bin. 
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8 X lO^'^ erg/s, which separates MIR-weak (small open symbols) and MIR-strong (filled symbols) sources 
(|Ogle et al. 2006). The large open symbols show the mean values of log-Rdr and Ln^MHz for the quasars 
and radio galaxies in four £178 mhz luminosity bins (see Table [5]). 



In order to quantify these trends, we divided the 
objects into four luminosity bins in ii78MHz- For 
each luminosity bin we then calculated the mean 
log i?dr and the mean Lusmhz for the radio galax- 
ies and the quasars separately. The results are in- 
dicated in Figure [S] and summarized in Table [51 
These numbers confirm that the ratio of i?dr be- 
tween radio galaxies and quasars changes with lu- 
minosity and that there seems to be a maximum 
difference at a luminosity ii78MHz of a few times 
lO'^^erg/s/Hz, where most sources belong to the 
intermediate redshift (0.5 < z < 1.0) sample. 
Because the quasars show a fairly continuous de- 
crease in i?dr over all Lusmuz (or z) covered here, 
the variation in the i?dr ratio is probably mostly 
due to the radio galaxies. 

One possible explanation of the trends we ob- 
serve is that the fraction of the total power emitted 
by the radio lobes increases with I/i78MHz- Be- 
cause for our flux-limited sample an increase in 
radio luminosity on average translates into an in- 
crease in redshift, one could imagine increased lobe 
power due to (on average) denser intergalactic en- 
vironments at higher redshifts. Such effects, how- 
ever, should apply to radio galaxies and quasars 



alike and cannot easily explain the i?dr upturn for 
the most luminous radio galaxies. 

The apparent _Rdr decline in Figure [8] could in 
principle be caused by the simple fact of plotting 
l/Li78MHz versus Li78MHz, but this again would 
not explain the rising i?dr for the most luminous 
radio galaxies. Therefore we assume that there is 
a real physical effect dependent on luminosity and 
differing for quasars and radio galaxies. 

One example of a possible physical effect is that 
the dust covering angle may be reduced at high 
luminosities, for example by destroying dust or 
pushing it farther away from the nucleus. This 
would-be consistent with the receding-torus model 
(e.g. Lawrence 1991). A decline of Rdi- with lumi- 
nosity is expected in such a model if the height of 
the torus remains constant. However, the height 
could be blown up by, for instance, accompany- 
ing starbursts, a nd this could also be correlated 
with luminosity ([Haas et al.ll2004) . But even in 
this picture the behaviour in i?dr for the most lu- 
minous radio galaxies remains a puzzle. 

An alternative explanation could involve the ef- 
fects of additional absorption of the radio galax- 
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Table 2: Average _Rdr for quasars and radio galaxies in four luminosity bins 



sample 
(ii78MHz in ergs"^ Hz" 


') 


number of 
sources 


(l0gi?dr) 


{fl.di>QSRs 
{fldr)RGs 


il78MHz < 1 X 10-^5 




QSR: 7 
RG: 11 


2.41 ±0.16 
2.03 ±0.18 


2.4 ±1.3 


1 X 10^5 <Li78MHz <5X 


1035 


QSR: 10 
RG: 18 


1.92 ±0.09 
1.10 ±0.09 


6.6 ±1.9 


5 X 1035 < Li78MHz < 1 X 


1036 


QSR: 11 
RG: 10 


1.67 ±0.07 
1.24 ±0.15 


2.7±1.0 


^178 MHz > 1 X 1036 




QSR: 8 
RG: 8 


1.53 ±0.07 
1.31 ±0.07 


1.7 ±0.4 



Note. — The average R^r is taken as (i?dr) = 10^'°^ ^■J'' with R^,- = i^Li^,i5;jm/^'ii/,i78MHz- MIR-weak sources with 
J^^i^jlS/jm = 8 X 10'*'^ erg/s have been excluded. Including the MIR-weak radio galaxies into the statistics does not change 
the trends; the (logiJdr) changes to 1.40 ± 0.16 and 0.94 ± 0.12 for the two lower luminosity bins, respectively. There are no 
MIR-weak objects in our sample for the two higher luminosity bins. 



ies' 15 /im continuum emission by dust in the host 
galaxy, but outside the torus. A high dust column 
density in the host would be required for such be- 
havior. In a situation where the dust (including 
the torus) is pushed to larger radii due to the in- 
creased luminosity of the nucleus, the line-of-sight 
dust column density in the host might become too 
small to make significant contributions to the to- 
tal absorption of the continuum emission. Mid- 
and far-IR observations of the optically-selected, 
mostly radio quiet Palomar-Green (PG) quasars 
indicate that the cold end of the dust tempera- 
ture distribution of the most luminous PG quasars 
at 1 < z < 2 is surprisingly warm compared 
with P G quasars of low and intermediate lumi- 
nosity (JHaas at al.ll2003[ ). This suggests that the 
amount of colder dust in these (hyper)-luminous 
IR sources is relatively small and that the AGN 
effectively heats much of the host galaxy's dust. 
If this scenario applies for our high-z 3CRR radio 
galaxies as well, then the total line-of-sight dust 
column density may be too small to be able to ab- 
sorb as much as one would expect from a scaled- 
up version of lower-luminosity AGNs, which are 
located in hosts relatively richer in dust. In other 
words, if the AGN power relative to the host's 
dust mass exceeds a critical threshold, the maxi- 
mum possible extinction declines, hence the ratio 
of i?dr(QSR)/i?dr(RG) declines as weh. 



Most likely there will not be a single explana- 
tion for the observed trends but a combination of 
effects will have to be considered. Additional in- 
sights may come from extending the luminosity 
range by including Herschel observations of the 
highest redshift 3CRR sources, or by including 
further constraints from a Chandra X-ray study 
of the present sample (B. Wilkes, in preparation). 

3.4. Emission Lines 

Four prominent emission lines can be identified 
in the average spectra (Figure [2|): the low ion- 
ization [Ne II] line, the medium ionization species 
[S IV] and [Ne III], and the high ionization AGN 
tracer [Ne V]. We determined emission-line fluxes 
only from the average spectra, but their uncertain- 
ties still are large. This is for two reasons. The 
individual spectra in which emission lines can be 
identified show a wide range of emission-line equiv- 
alent widths for each type of object. This disper- 
sion among individual objects leads to consider- 
able uncertainty for the line fiuxes in the average 
spectra. In addition, the measured emission-line 
fluxes in the average spectra depend strongly on 
the placement of the continuum. Especially for 
weak or noisy features on top of a complex con- 
tinuum, this results in large errors. For these rea- 
sons, conclusions from the emission lines should be 
considered tentative. With these caveats in mind, 
within the uncertainties the emission line ratios 
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Table 3: Emission-line lum.inosities. 



A 


[SIV] 
10.5 /zm 


[Ne II] 
12.8 ^m 


[NeV] 
14.3 Aim 


[Ne III] 
15.5 Aim 


QSRs 
RGs 


2.5 
2.2 


1.2 
1.5 


1.9 
1.6 


4.6 
4.2 













Note. — The emission-line luminosities are measured in 
the average spectra and given in 10^^ erg/s. The uncertainties 
are ~25% for the brightest line [Ne III] and ^50% for the other 
three species (see text). 



for radio galaxies and quasars are similar. More- 
over the line luminosities are comparable as well 
(Table [S]). 

The spectra cover the position of the 11.3 /zm 
PAH feature, but this star formation tracer was 
not detected in the average spectra. This suggests 
that the MIR emission in our objects is strongly 
dominated by the AGN and that any contribu- 
tions from star formation are small compared to 
AGN emission. 3CR radio sources at z < 1 also 
have less star formation than optically or NIR se- 
lected radio-qui et QSOs at comparable redshifts 
(|Shi et al.l 120071 ). In fact, in most lower-redshift 
SCR sources, PAH features remain undetected. 
This is consistent with the finding that many of 
the radio sources are hosted by rather q uiescent 
elliptical galaxies fe.g.. lFlovd et al.ll2008l ). 



4. Conclusions 

The 3CRR sample at 1 < z < 1.4 represents 
the most radio-luminous steep-spectrum quasars 
and double-lobed radio galaxies for which the rest- 
frame MIR 9-16 /im wavelength range is accessible 
with the Spitzer Infrared Spectrograph. The sam- 
ple is unbiased with respect to orientation. 

1) The quasar SEDs are very similar to each 
other in shape. They show a strong NIR 
and MIR continuum, silicates in emission, 
and an emission bump at 2-5 /ini. 

2) The mean radio galaxy spectrum (8 < 
Arcst < 16 /im) is consistent with a quasar 
spectrum reddened by a dust screen of 
Av ^ 20 mag both in spectral shape and 
in flux level. 



3) Radio galaxy extinctions indicated by their 
2-5 /im SEDs are higher than the extinction 
derived from the silicate absorption and the 
(rest) 15 /im flux. This can be explained 
if the component responsible for the NIR 
bump in the quasar SED suffers more ex- 
tinction than the MIR emission. 

4) All sources in our sample with silicate ab- 
sorption have small core fraction (<10~*), 
while objects with silicate emission often 
show stronger cores. This strongly supports 
the idea that the observed extinction in our 
sample is related to the orientation of the 
source. 

5) The average spectra show similar emission- 
line ratios for radio galaxies and quasars, 
and no PAH features are detected. This ar- 
gues in favor of similar intrinsic AGN power 
and star-forming activity in quasars and ra- 
dio galaxies, as predicted by unified schemes. 

6) The comparison of our high-z sample with 
lower redshift 3CRR sources reveals a steady 
decline of the quasar i?dr with increasing lu- 
minosity, and a partial decline, followed by a 
high-luminosity upturn, of the radio galax- 
ies' i?dr- Possible explanations for these 
trends may involve an increased fraction of 
power emitted by the radio lobes due to 
different galactic environments at higher z 
and/or different dust geometries due to the 
higher nuclear luminosities (e.g distance of 
the dust from the nucleus, solid angle of the 
dust exposed to the radiation, or total dust 
column density on the line of sight). 

While the MIR spectra provide valuable insight 
into the nature of the MIR emission and the rela- 
tion between the most powerful type-1 and type- 
2 radio AGN, we expect further constraints from 
Chandra X-ray and Herschel far infrared/sub-mm 
observations of this sample. 
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